Distribution of reflection coefficients in absorbing chaotic microwave cavities 
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The distribution of reflection coefficients P(R) for chaotic microwave cavities with time-reversal 
symmetry is investigated in different absorption and antenna coupling regimes. For all regimes the 
agreement between experimental distributions and random-matrix theory predictions is very good, 
provided both the antenna coupling T a and the wall absorption strength T w are taken into account 
in an appropriate way. These parameters are determined by independent experimental quantities. 
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Wave scattering by chaotic and weak disordered sys- 
tems has motivated a rather intensive research activ- 
ity. The subject is common to several areas of physics, 
ranging from nuclear, atomic, molecular, and mesoscopic 
physics to classical wave scattering, like microwaves, 
sound, and light (see Refs. 0,0] for a review). The 
striking feature shared by all such systems, provided they 
have a chaotic underlying classical dynamics, is that they 
show universal transmission fluctuations. These fluctua- 
tions are successfully described by random-matrix theory 
0, 0] . Very recently a comprehensive treatment of ab- 
sorption, ubiquitous in experiments, was developed for 
systems without time-reversal symmetry (TRI) 0, 
For systems with TRI, however, there is a rigorous theory 
only in limiting cases Q. 

An analytical expression for the distribution P(R) of 
the reflection coefficient R in the presence of absorption 
was derived for an arbitrarynumber of open channels in 
the weak absorption limit [7j. In the strongabsorption 
limit P(R) reduces to a simple exponential [8j. In these 
two works perfect coupling between the channels and the 
scattering region is assumed. Furthermore, non-resonant 
backscattering processes U fiol Hlj , which are present in 
most experimental situations, have been until now often 
neglected in cases that absorption is taken into account. 
It is noteworth y th at only few of these works are experi- 
mental [II El III . 

The purpose of this letter is to present experimental 
evidence that random-matrix theory provides a quanti- 
tative understanding of the universal reflection fluctua- 
tions in chaotic systems. More specifically, we show that 
this can only be achieved if both the coupling, so far 
mostly overlooked, and absorption strengths are properly 
taken into account. To this end we measure P(R) in mi- 
crowave cavities from weak to strong absorption regimes. 
A quantitative agreement between experiment and the- 
ory is observed for all the cases we have studied. This is 
remarkable since all theoretical parameters are directly 
obtained from averaged experimental quantities and not 
from a fit of P(R) to the data. 

To study the dependence of P(R) on absorption and 
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FIG. 1: Sketch of the microwave billiards, (left) Half Sinai 
billiard (ai = 43 cm, bi = 23.7 cm, height hi = 7.8 mm). The 
half disk (d = 12 cm) was moved along oi. For the second 
set of measurements one wall was coated with an absorber, 
(right) Disordered billiard (d2 = 34 cm, 62 = 24 cm, height 
h,2 = 8 mm) [Ig. The billiards are drawn at scale. 



coupling strengths, we used three different flat mi- 
crowaves cavities and measured the reflection. For the 
half Sinai billiard, Fig. QJleft), 57 reflection spectra were 
measured by sliding the half circular inset along the wall 
over 28.5 cm in steps of 0.5 cm. This is a practical way 
to improve the statistics, notwithstanding the correla- 
tions among spectra. For all cases only one antenna is 
coupled to the cavity. For the second measurement, one 
wall of the half Sinai billiard was coated with an absorber 
|l5|. The third measurement was performed on a disor- 
dered billiard [16| , where 20 brass cylinders of 5 different 
sizes have been located at random within the billiard 
(Fig. fright)). The data were collected on a 5mm grid 
with 696 antenna positions. We do not consider antenna 
positions closer than 15 mm to the billiard walls, so only 
533 measurements are included in the statistical analysis. 
The reflection coefficients R were measured with vector 
network analyzers. Typical absorption spectra for the 
three cases are shown in Fig. [21 

Figure 13 shows the mean reflection coefficient (R) as a 
function of the frequency v for the different cavities. (• • •) 
denotes averaging over a frequency range Ais (see Fig. 
and over different cavity geometries (half Sinai billiard) 
or different antenna positions (disordered billiard). Five 
representative windows, where (R) is approximately con- 
stant, are chosen for a detailed statistical investigation. 
The windows correspond to the shaded areas in Fig. [21 
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FIG. 2: Typical absorption spectra for different absorption 
regimes: (a) Weak absorption and weak coupling - disordered 
billiard, (b) intermediate absorption and coupling - half Sinai 
billiard, and (c) strong absorption and perfect coupling - half 
Sinai billiard with absorber. 



We shall now discuss how P(R) is obtained using 
random-matrix theory The scattering matrix S is 
unitary and has a dimension N, the number of open 
channels. Since absorption opens new channels, the ex- 
perimentally accessible 7V cxp -channel sector of 5, which 
we call S 1 , becomes subunitary. As a result, the eigenval- 
ues Ri, i?2, ■ ■ ■ , Rn cxp of SS^ lie between and 1 p]. In 
the whole frequency range the diameter of the antenna 
is small compared to the wavelength implying a single 
coupled channel, i. e., there is a single measurable reflec- 
tion coefficient R. The resonant S matrix for quantum 
systems |l7j can be used to describe cavity microwave 
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FIG. 3: Average reflection (R) as a function of the frequency 
v for: half Sinai (triangles), half Sinai with absorber (squares) 
and disordered (diamonds), cavities. The shaded boxes mark 
the cases subjected to a detailed analysis. The frequency 
average is A^=450 MHz. 

scattering 0, 0] , and reads 

S{E) = 1 - 2itiW\E-H + inWW^)~ 1 W . (1) 

The quantum mechanical energy E is related to the elec- 
tromagnetic frequency v by Weyl's formula. In the high- 
frequency limit E ~ v 2 , whereas at low frequencies the 
dispersion relation depends on the cavity shape. H de- 
scribes the cavity quasi-bound spectrum, and W contains 
the coupling matrix elements between the resonances and 
the scattering channels. These are the absorption chan- 
nels in the cavity N w 3> 1 and the microwave antenna 
N a = 1 ■ These absorption channels act as additional 
"fictitious" waveguides attached to the cavity |2cj . 

In our statistical treatment we take the matrix H as 
being a member of the Gaussian orthogonal ensemble, as 
usual for time-reversal invariant chaotic systems. Since 
the H matrix is statistically invariant under orthogonal 
transformations, the statistical properties of 5* depend 
only on the mean resonance spacing A (determined by 
H) and on the traces of W^W. The average scatter- 
ing properties are best characterized by the transmission 
coefficients of the antenna t a and of the fictitious waveg- 
uides t w . These coefficients are defined as t c = I — IS^d 2 , 
where c is an arbitrary channel and 777 stands for an 
ensemble average. The ergodic hypothesis identifies the 
later with the experimental frequency averages. Weak 
coupling corresponds to t c — > 0, or direct reflection from 
channel c (no flux entering the cavity) , whereas t c — > 1 is 
the limit of perfect coupling. Hence, our model has two 
parameters, namely, the antenna coupling T a = t a and 
the absorption strength T w — N w t w . 

Analytical results for P(R) in the case of time-reversal 
invariant systems are known only for the limiting cases 
of very weak (T w -C 1) or strong absorption (T w ^> 1) 
with perfect antenna coupling (T a = 1). For most cases 
of interest one has to rely on numerical simulations of 
the resonant S matrix as given by Eq. (Q. This is done 
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billiard 


A^/GHz 


(R) 






disordered 


3.5- 4.0 


0.981 


0.011 


0.49 


half Sinai 


4.0- 5.0 


0.906 


0.116 


0.56 


half Sinai 


8.0- 9.0 


0.517 


0.754 


2.42 


half Sinai 


14.7-15.7 


0.185 


0.989 


8.40 


absorbing wall 


13.5-15.5 


0.039 


0.998 


48.00 



TABLE I: Mean reflection coefficient (R), antenna coupling 
T a , and wall absorption T w for the cases marked in Fig. [3] 
Av is the frequency range. 

as follows: first, T a and T w determine W 0]; second, 
by picking H from the Gaussian orthogonal ensemble we 
compute S. In this study we generate 10 5 realizations of 
H with dimension M = 201 and fix N w — 200 for each 
case of interest. These values are justified as follows: The 
exact value of N w is not of relevance, provided N w ^> 
1, as demonstrated in Ref. '21\. The statistical theory 
requires M 3> 1, which is the case in our experiment. 
Finally, for technical reasons M > N w [l7| . 

The transmission coefficients T a and T w are obtained 
directly from the experimental data as follows: First, we 
use t a = 1 — K^a)! 2 to determine T a directly from the 
measured values of S aa {y)- T w is then fixed by adjust- 
ing the theoretical mean reflection coefficient (R) to the 
experimental one. In the low frequency regime, the val- 
ues of (Saa(v)) and (R) cannot be extracted correctly 
from the experimental data since the background of the 
reflection spectrum varies slowly with v and deviates 
from unity, as shown in Fig. |2fa)- These variations of 
the background cause a broadening of the distribution 
and a shift of its maximum. This problem can be over- 
come as follows: Since in this regime the resonances are 
nearly isolated, each of them can be fitted by a complex 
Lorentzian. Using these fits a corrected spectrum can 
be generated with a perfect background and suppression 
of noise. From this spectrum (S aa ), (R), and P(R) are 
obtained. The values for the mean reflection coefficient 
(R), the antenna transmission T a , and the wall absorp- 
tion T w for the five regions under investigation are shown 
in Table IJ 

We first discuss the regime of strong absorption. The 
theoretical distribution is given by || 

P(R) = (Ry 1 e- R/{R \ (2) 

which no longer depends explicitly on the number of 
open channels. Experimentally, the system with the 
largest absorption is the half Sinai billiard with absorber 
as seen in Fig. |3J A typical spectrum is presented in 
Fig. Gfc) . The corresponding experimental distribution 
P(R) is shown in Fig. together with the prediction of 
Eq. (J2J and the numerical simulation. The theory con- 
tains no free parameters since the only inputs, namely T a 
and (R), were independently extracted from the experi- 
mental data. A strong variation of T w in the investigated 
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FIG. 4: Distribution of the reflection coefficient R for the half 
Sinai billiard in the strong absorption limit (13.5-15.5 GHz). 
The solid line corresponds to Eq. and the dotted one to 
the numerical simulation. 
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FIG. 5: Distribution of the reflection coefficients, P(R), for 
weak coupling strength and small absorption (3.5-4.0 GHz), 
obtained from the Lorentzian fit corrected spectrum (for de- 
tails see text). For comparison the distribution obtained from 
the uncorrected spectrum is shown in light grey. The inset en- 
larges the region of R between 0.95 and 1.005. The dotted 
line is obtained from the random-matrix simulation. 

frequency range causes the discrepancy between the ex- 
periment and the simulation. 

Let us now turn to the weak absorption regime, T w < 
1, at low frequencies. Unfortunately, the available ana- 
lytical result |3 cannot be used, since it assumes perfect 
coupling, whereas in the present experiment, the coupling 
is weak at low frequencies (see Table [IJ . A similar distri- 
bution for nearly perfect coupling has been presented in 
Ref. [IJ, but no detailed analysis is performed. In Fig. [S] 
the distribution of P(R) is shown for the disordered cav- 
ity in the frequency region of 3.5 to 4.0 GHz, obtained 
from the Lorentzian-fit corrected-spectra. This is be- 
yond the localization-delocalization transition observed 
at about 3 GHz for the disordered cavity of Fig. right) 
(a selection of eigenfunctions can be found in Ref. |lq). 
The results from the numerical simulation arc plotted for 
comparison with almost perfect agreement. 

The experimental distributions P(R) for intermediate 
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FIG. 6: Distribution of the reflection coefficients P(R) for 
the intermediate absorption regime (half Sinai billiard). The 
data were collected at different frequency windows marked in 
Fig. 03 The dashed lines stand for the random-matrix simula- 
tions. In (a) the experimental distribution was obtained form 
the Lorentzian fit corrected spectrum. 



absorption and coupling regimes are presented in Fig. [5] 
The data are taken from the half Sinai billiard depicted 
in Fig. [IJleft). A typical spectrum for the intermedi- 
ate regime is shown in Fig. [2b). In all cases, a very 
good agreement between the experimental and theoreti- 
cal P(R) is found. 

It was essential to extract the coupling and the absorp- 
tion strength from two independent experimental quan- 
tities. The coupling strength depends on the phase vari- 
ations of S through (S), whereas the absorption strength 
depends on the average (R) . We obtain a good agreement 
between experiment and theory, comparable to the one 
displayed in Fig-EI by assuming perfect coupling (T a = 1) 
and thus overestimating the absorption. Such agreement 
can lead to a wrong understanding of the actual physical 
process at hand. This is an important message to bear 
in mind when analyzing quantum mechanical processes, 
where in most cases the phase of S cannot be measured. 



We have shown that microwave experiments with 
chaotic and disordered billiards are excellent tools to in- 
vestigate the universal distribution of reflection coeffi- 
cients with absorption. We studied a broad variety of 
regimes, ranging from weak coupling and weak absorp- 
tion to perfect coupling and strong absorption. In all 
cases a very good agreement between experiment and 
theory was found. This is particularly striking, since the 
theoretical input parameters were fixed independently by 
the experiment. 
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